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Abstract 

A kinetic equation describing the nonlinear evolution of intense electromagnetic pulses in 
electron-positron (e-p) plasmas is presented. The modulational instability is analyzed for a rel- 
ativistically intense partially coherent pulse, and it is found that the modulational instability is 
inhibited by the spectral pulse broadening. A numerical study for the one-dimensional kinetic 
photon equation is presented. Computer simulations reveal a Fermi-Pasta-Ulam-like recurrence 
phenomena for localized broadband pulses. The results should be of importance in understanding 
the nonlinear propagation of broadband intense electromagnetic pulses in e-p plasmas in laser- 
plasma systems as well as in astrophysical plasma settings. 

PACS numbers: 52.38.-r, 52.27.Ny, 52.65.Rr 
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I. INTRODUCTION 



Electron-positron plasmas are believed to be an important ingredient of the early universe 
and in astrophysical objects such as pulsars, supernova remnants and active galactic nuclei, 
and in gamma-ray bursts jjij. In such extreme environments, the electron-positron pairs 
may be created by collisions between particles that are accelerated by electromagnetic and 
electrostatic waves and/or by gravitational forces. In pulsar environments, there is also a 
possibility of pair creation via high-energy curvature radiationphotons that are triggered by 
charged particles streaming along the curved magnetic field |2j, with a resulting collection 
of positrons at the polar caps of the pulsar Q, Q]. High-energy laser-plasma interactions 
and fusion devices on Earth also constitute a source of electron-positron plasmas. Experi- 
ments with petawatt lasers (with intensities exceeding 10 20 W/cm 2 ) have demonstrated the 
production of MeV electrons and evidence of positron production via electron collisions 
Positrons are also believed to be created in post-disruption plasmas in large tokamaks 
through collisions between MeV electrons and thermal particles jzj. 

Various collective phenomena in electron-positron plasmas have been demonstrated in the 

n 

laboratory, such as the creation of wake fields by ultra-relativistic positron beams |8| and 
two-stream instabilities in streaming electron-positron plasmas [9]. The natural saturation 
mechanism for a two-stream instability is particle trapping in which electrons and positrons 
are trapped in the electrostatic potential of the large-amplitude wave. Theoretical and 
numerical studies have shown that solitary electron/positron holes can exist both in a non- 



relativistic pair plasma [10|, [111, [12j and in a relativistically hot electron-positron-ion plasma 
[13I 1 . Intense electromagnetic radiation introduces a new nonlinearity in the electron-positron 
plasma via the relativistic mass increase of the plasma particles in the quivering electric field 
of the electromagnetic waves Q, Q] . It has been demonstrated theoretically that relativistic 
electromagnetic solitary waves can exist in a cold electron-positron plasma, if there is an 
inclusion of massive ions [3]. Several authors 12] have also considered the nonlinear 
propagation of small amplitude electromagnetic waves and Alfven vortices in a strongly 
magnetized electron-positron plasma. 

Over recent years, the laser pulse intensity has been strongly increasing, and over the 
next five to ten years it is expected to reach focal intensities of ~ 10 26 W/cm 2 [19J. Thus, 
the interaction between plasmas and intense laser pulses is currently of great interest [20, 
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24j . As laser intensities keep increasing, laboratory plasma experiments will 



experience the effects of the nonlinear quantum vacuum [25[, such as pair creation and 
photon splitting. Theoretical studies of nonlinear laser-plasma interactions mostly make 
the simplified assumption of a coherent, monochromatic wave propagation in plasmas. It is 
then relatively easy to derive analytic expressions for various instabilities (e.g. relativistic 
Raman forward and backward scattering instabilities and the modulational instability in an 
electron- ion plasma [2^) and the formation of envelope electromagnetic solitons 2?]]. 

However, the assumption about a monochromatic wave is not always true, e.g. if the 
electromagnetic waves are excited via turbulent processes. Then, we would have a distri- 
bution of waves with different wavelengths and frequencies that are only partially coherent. 
Theories for photon acceleration of partially coherent waves have been developed in the 

n n 

framework of nonlinear optics [28[ and are based on the Wigner's kinetic description [29]. 
The latter has been used to investigate the nonlinear instability of random phased small 
amplitude electromagnetic waves in a nonlinear dispersive medium 3fl| . 

In this paper, we consider the nonlinear propagation of relativistically intense broadband 
laser pulses in an electron-positron plasma. In order to analyze the properties of such 
partially coherent laser pulses, we perform a Wigner analysis of the modified nonlinear 
Schrodinger equation, and obtain a wave kinetic equation for the photon quasi-particles. 
The resulting kinetic equation is investigated both analytically and numerically. Specifically, 
we present the modulational instability and the dynamics of partially coherent intense laser 
pulses. It is found that the partial coherence of the laser pulse yields a reduced growth rate 
for the modulational instability. The fully nonlinear evolution of the wave kinetic equation 
is investigated by means of numerical simulations. The latter show a complex dynamics of 
broadband laser pulses in that the laser envelope contracts and disperses in a quasi-periodic 
manner, somewhat similar to the Fermi-Pasta-Ulam recurrence phenomenon encountered for 
the nonlinear Schrodinger equation. The present results should be useful for understanding 
the nonlinear propagation of broadband intense electromagnetic pulses through pair plasmas 
such as those in laser-plasma systems and in astrophysical settings. 
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II. GOVERNING EQUATIONS AND THE MODULATIONAL INSTABILITY 



Let us consider the nonlinear propagation of intense circularly polarized electromagnetic 
waves in pair plasmas without ions. Accounting for the relativistic mass increase of the pairs 
and quasi-stationary density fluctuations driven by the relativistic ponderomotive force, 
I. derived the fully nonlinear equation 
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l exp[,3(l- V l + |a| 2 )] 



\a 



2 



0, (1) 



which shows the evolution of the normalized vector potential in pair plasmas. Here a = 
eA/mc 2 is the normalized vector potential, A is the amplitude of the circularly polarized 
vector potential, m is the electron mass, e is the magnitude of the electron charge, c is the 
speed of light in vacuum, (3 = 2f3 e P P / (Pe + P P ) is the dimensionless temperature parameter 
with (3e iP = (c/vTe, P ) 2 , VTe, P = {T^p/m) 1 ! 2 is the thermal speed, T e (T p ) is the electron 
(positron) temperature, we have normalized the time and co-moving spatial variables by 
ujq/lu 2 and c/lu p , respectively, lu p = (4irn e 2 /m) 1//2 is the electron plasma frequency, n is 
the unperturbed electron density, and ooo is the central wave frequency. Numerical analysis of 
Eq. (1) revealed that weakly modulated electromagnetic (em) pulses would undergo collapse, 
leading to strong intensification of the localized em pulses. 

Even though Eq. (PJ) contains a saturation nonlinearity, halting of collapse could be ob- 
tained by spectral broadening techniques, well-known in inertial confinement fusion. More- 
over, in many applications of such an equation, for example astrophysical systems, the pulses 
can be partially coherent. Thus, understanding the dynamics and stability of electromag- 
netic pulses taking incoherence effects into account may be important for practical purposes. 

A canonical way for analyzing effects of partial coherence is to use the so called Wigner 
function, defined as the Fourier transform of the two-point correlation function of the vector 
potential a according to [3, |3] 

p(t,r,k) = ^L J dee*-V(t,r + £/2)a(t l r-£/2)) l (2) 

where the angular bracket denotes the ensemble average. The Wigner function represents a 
generalized distribution function for quasi-particles, in this case photons. From the definition 
(J2J), one finds the relation 

J(t,r)= / dkp(t,r,k), (3) 



where I = (|a| 2 ). Applying the time derivative to the definition (J2J) and using Eq. (JTJ), we 
obtain the kinetic equation 



dp , „ 2exp[/3(l - Vl+T)] . A - - \ 

for the quasi-particles. Here the sin-operator is defined in terms of its Taylor expansion. 
Keeping only the first term in the latter, which corresponds to the long wavelength limit, 
we obtain a photon kinetic (or the Liouville) equation 



exp[/3(l- VTT7)] 



V feP = 0. (5) 



VT+i 

However, the photon kinetic limit is only valid for weak spatial variations in p, and produces 
a modulational instability growth rate which is unbounded as the wavenumber of the photons 
increases. Thus, for a broad spatial spectral distribution of photons the dynamics of the 
Liouville equation is highly unstable and strongly deviates from the full Wigner dynamics, 
and will therefore not be analyzed further here. 

Equations (JHJ) and (jlj) constitute a complete description of partially coherent nonlinear 
photons in an electron-positron plasma. 

Next, we perform a perturbation analysis of Eq. (jlj). Letting p = po(k) + piexp(iK • 
r — iflt), where \pi\ <C p , we linearize Eqs. © and (jlj) with respect to the perturbation 
variables. From Eq. (J3J) we then obtain 

-iQ Pl + ^k • K Pl - sin ■ V fc ) p = 0, (6) 

where 



exp[/?(l- yO+To)] 

Combining (5) with Eq. (JHJ) we obtain the nonlinear dispersion relation 

dU(Io) f dh Po(k-K/2)-p (k + K/2) (g) 



dJ J Q-kK 

which is valid for partially coherent intense laser pulses in an electron-positron plasma. 

In the case of a monochromatic pulse, the background distribution function satisfies 
Po(k) = Iq5(\l — ko) for some wavevector ko. The dispersion relation (jHJ) then reads 

1 . , T dU{IJ 1 ' 



tt = K • k ± 



4 dl 



(9) 
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which is agreement with the results found in Ref. |31|. The growth rate T 
by 

r 



2 V 1 + 4 



lV2 



K 2 - -K 4 
4 



-iQ is given 



(10) 



In order to simplify the perturbation analysis, we introduce the assumption of one- 
dimensional partially coherent photon propagation along the z-axis. We investigate the 
case of partial coherence using the background Lorentz distribution 

I A 



Po(fc) 



(11) 



7T (k - k ) 2 + A 2 ' 

where k = k z . The distribution (jllj) corresponds to a partially coherent phase of the vector 
potential do, giving rise to a spectral broadening of po with a width A. From the dispersion 
relation (El) we then obtain 



-iKA + Kk n ± 



4 



r dU(Io) k2 
1q ~Tt — K 



1/2 



(12) 



thus giving Q in the limit A — ► 0. We see that the effect of the spectral broadening is to 
reduce the growth rate according to 



-KA 



2 [P + 



U{h) 
l + /o 



K 2 



1/2 



(13) 



We note that as A — > 0, we retrieve the expression (fTU|) . 

A comparison between the coherent and incoherent modulational instability growth rates 
is presented in Fig. 1. We note from Eq. (|13|) that the growth rate is larger than zero only 
if A < (/o/2) 1/2 [/9 + U(I )/(1 + / )] 1/2 , and in the small amplitude limit J « 1 we have 
the condition A < (IqP/2) 1 ^ 2 for the modulational instability The Liouville equation is 
obtained from the Wigner's kinetic equation when the product of the spatial length scale L 
and the spectral width A is large enough. For ultra-short laser pulses with a focal intensities 



of J 



Jq\A\ 2 /c ~ 10 21 W/cm 2 , wavelength Ao ~ 800 nm, and spectral width 5X ~ 50 nm 
(relevant for e.g. the Astra laser and its upgrade Astra Gemini at the Rutherford Appleton 
Laboratory, Oxfordshire, U.K.), we have the coherence length l c = Ag/5A ~ 13 /im. Since 
A ~ c/£ c u p we have the modulational instability criteria (% / M C rit)(c/t>Te) 2 > m 2 c 5 /47re 2 X£ 2 , 
where n crit = mio\j 'Ane 2 . For the above intensity the criteria for the modulational instability 
to take place then becomes 2.1 x 10~ 6 < {n /n CT [ t )(c/vT e ) 2 , and if the electron temperature 
is T e ~ 5keV, we find that n > 4 x 10~ 22 n crit ~ 10 -2 cm -2 . 
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III. NONLINEAR DYNAMICS OF BROADBAND INTENSE ELECTROMAG- 
NETIC PULSES 



In order to understand the long term behavior of modulationally unstable broadband 
intense electromagnetic pulses in pair plasmas, we carry out the numerical analysis of the 
kinetic photon equation (@J in the one-dimensional case. Accordingly, we numerically solve 
the set of equations 



dt x dx 



2exp[/3(l - VT+J)] .(id d 

2 dx OK 



p 



o. 



and 



I(t,x) = 

For the numerical solutions of Eqs. 



dk x p(t, x, k x 



(14) 



(15) 



Vlasov-Poisson system of equations 



and (fT3j) . we use a Fourier method for the 



32j-L!i|- which is slightly modified to solve the Wigner 



equation. For this purpose, we use the Fourier transform pair 

/oo 
drjp(t,x,ri)e- ivk *, 
-oo 

1 f°° 

p(t,x,7]) = — dk x p(t,x,k x )e ir]kx 
to obtain the Fourier transformed version of Eqs. and (fTK|) . 



dp . d 2 p ^ . 
dt dr/dx ^ 



it] d 
~2dx 



2exp[/?(l- VT + 7)] 



and 



(16) 
(17) 

(18) 
(19) 



I(t,x) = 27rp(t J x,r]) v=0} 

respectively. The Fourier transformed system of Eqs. (fTHj) and ()19|). which is equivalent 
to Eqs. (fTl| and (fTB^) . is solved numerically in a periodic box in x space, where the x 
derivatives are approximated with a pseudo-spectral method. In this method, the differential 
operator d/dx turns into a multiplication by iKj and the operator sm[(ii]/2)d/dx} turns 
into a multiplication by sin(— rjKj/2), where Kj = 2irj/L is the spatial wavenumber, j = 
0, ±1, ±2 . . . and L is the length of the spatial domain. We note that the integral over k x 
in Eq. (fTB^) is transformed into a simple evaluation of p at rj = in Eq. ()19|) . As an initial 
condition for the simulation, we take the Fourier transformed Lorentz distribution (|TT| 

h 



Po(v) 



2tt 



exp(— A|?y| + ikor}), 



(20) 



and the same parameters as in Fig. 1. In Figs. 2 and 3, we use Jo = 1/4, A = 0.1, and 
ko — 0.1 as initial conditions. A small amplitude noise (random numbers of order 10~ 3 ) is 
added to the initial condition to give a seed for the modulational instability. 

We display the temporal evolution of the em wave intensity / in Fig. 2. Here we see 
an initially linear growth phase and a wave collapse that takes place around t = 90. The 
wavelength of the collapsing pulse is approximately A = 10, and thus the most unstable mode 
is k* = 27r/10 w 0.6, in good agreement with the theoretical prediction in Fig. 1. After the 
initial pulse collapse, the solution shows a chaotic behavior, where some localized envelopes 
show a semi-periodic behavior, and we also see a sequence of merging and splitting of wave 

am recurrence phenomenon 



and has been explained 

n n 



groups. This behavior is somewhat similar to the Fermi-Pasta-U 
|34|, which has been observed in nonlinear optics experiments 

theoretically in the framework of the nonlinear Schrodinger equation |36j, |37|. In Fig. 3, 
we show Wigner's distribution function p at different times. (It is obtained by numerically 
inverse Fourier transforming p to obtain the real-valued p.) We note that the distribution 
function shows a complex dynamics and takes both positive and negative values, in contrast 
to solutions of the Vlasov equation for particles that takes only positive values. In Figs. 4 and 
5, we repeated the simulation with the larger spectral width A = 0.3 in the initial conditions. 
We used the same intensity Iq = 1/4, as in the first simulation, and we used ko = 0.05. For 
this case, the temporal evolution of the pulse intensity, as shown in Fig. 4, exhibits a slower 
growth rate (compared to Fig. 2 with the smaller A = 0.1), and nonlinearly collapsing 
wave envelopes have smaller intensity maxima. The snap shots of the Wigner distribution 
function, depicted in Fig. 5, shows the initial instability and merging of nonlinear structures 
corresponding to merging wave envelopes in Fig. 4. Thus, a larger spectral width A leads 
to a slower and less violent dynamics of the wave envelopes that also have larger scale sizes 
than for the smaller values of A. We recall that the theoretical treatment predicts that there 
is a largest A above which the modulational instability vanishes, and we would instead have 
damping of waves similar as Landau damping for the Vlasov equation. 

IV. SUMMARY 

In summary, we have investigated the modulational instability and the nonlinear dynam- 
ics of partially coherent intense electromagnetic (em) waves in an electron-positron plasma, 
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taking into account the relativistic mass increase of the pairs as well as large scale density 
fluctuations that are created by the relativistic ponderomotive force of em waves. The dy- 
namics of broadband intense em waves is governed by a photon kinetic (or the Liouville) 
equation. The latter is analyzed to obtain a nonlinear dispersion relation, which admits the 
modulational instability growth rate for a Lorentzian distribution of partially coherent in- 
tense em waves. It is found that a broadband of em waves leads to a reduction of the growth 
rate. Furthermore, the photon kinetic equation has been numerically solved to understand 
the nonlinear dynamics of modulationally unstable broadband intense em pulses. Our sim- 
ulation results reveal the formation of localized em wave packets that show some similarity 
to the Fermi-Pasta-Ulam recurrence phenomenon. For larger values of the spectral width, 
the system shows a dynamics on a slower timescale and on larger length scales, in agreement 
with the linear theory. The present results should help to understand the nonlinear propa- 
gation of broadband intense electromagnetic waves in pair plasmas that appear in inertial 
confinement fusion schemes as well as in pulsar magnetosphere and supernovae remnants. 

This work was partially supported by the Swedish Research Council and the Deutsche 
Forschungsgemeinschaft . 
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Figure 1: The modulational instability growth rates plotted as a function of K, as given 
by the expression (JT3|) . We have used Iq = 1/4 and [3=1. The thick curve represents the 
coherent case with A = 0, the thin curve has a finite spectral width A of 0.1. The decreased 
growth rate due to spectral broadening can clearly be seen. 

Figure 2 (Color online): The intensity / as a function of x and t. The parameters used in 
initial are condition 1 = 1$ = 1/4, A = 0.1 and k$ = 0.1. 

Figure 3 (Color inline): The Wigner distribution function p as a function of x and k x , at 
times t = 0, t = 85, t = 100 and t = 150 (upper to lower panels), with the correspond- 
ing intensity in Fig. 2. For the initial conditions we used I = Iq = 1/4, A = 0.1 and k Q = 0.1. 

Figure 4 (Color online): The intensity / as a function of x and t. The parameters used in 
the initial condition are J = Jo = 1/4, A = 0.3 and ko = 0.05. 

Figure 5 (Color online): The Wigner distribution function p as a function of x and k x , at 
times t = 0, t = 300, t = 400 and t = 500 (upper to lower panels), with the corresponding 
intensity in Fig. 4. For the initial conditions we used I = Iq = 1/4, A = 0.3 and ko = 0.05. 
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